Pregnancy is linked to faster epigenetic aging in young women.

A central prediction of evolutionary theory is that energy invested into reproduction comes at the expense of somatic maintenance and repair, accelerating biological aging. Supporting this prediction are findings that high fertility among women predicts shorter lifespan and poorer health later in life. However, biological aging is thought to begin before age-related health declines, limiting the applicability of morbidity and mortality for studying the aging process earlier in life. Here, we examine the relationship between reproductive history and biological aging in a sample of young (20 to 22yo) men and women from the Cebu Longitudinal Health and Nutrition Survey, located in the Philippines (n = 1,735). We quantify biological aging using six measures, collectively known as epigenetic clocks, reflecting various facets of cellular aging, health, and mortality risk. In a subset of women, we test whether longitudinal changes in gravidity between young and early-middle adulthood (25 to 31yo) are associated with changes in epigenetic aging during that time. Cross-sectionally, gravidity was associated with all six measures of accelerated epigenetic aging in women (n = 825). Furthermore, longitudinal increases in gravidity were linked to accelerated epigenetic aging in two epigenetic clocks (n = 331). In contrast, the number of pregnancies a man reported fathering was not associated with epigenetic aging among same-aged cohort men (n = 910). These effects were robust to socioecological, environmental, and immunological factors, consistent with the hypothesis that pregnancy accelerates biological aging and that these effects can be detected in young women in a high-fertility context.

DNA methylation-based estimators of telomere length show low correspondence with paternal age at conception and other measures of external validity of telomere length.

In humans, DNA methylation (DNAm) based estimators of telomere length (TL) have been shown to better predict TL-associated variables (e.g., age, sex, and mortality) than TL itself. The biological significance of DNAm-based estimators of TL (DNAmTL) is unclear. In vitro DNAmTL shortens with cell replications, even when telomerase is maintaining TL. Telomerase is typically suppressed in humans, except in testes. Accordingly, sperm TL increases with age, and offspring with greater paternal age at conception (PAC) have longer TL. Thus, we expect that PAC associations with DNAmTL can shed light on whether in vivo cell replications in the presence of high telomerase activity (production of sperm) shorten DNAmTL or if PAC-lengthened TL causes lengthened DNAmTL. In a pre-registered analysis, using data from 1733 blood samples from the Philippines, we examined the association between paternal age at conception (PAC) and offspring DNAmTL. We did not find an association between PAC and DNAmTL but found a positive association of paternal grandfather's age at father's conception predicting grandchild's DNAmTL. In post hoc analyses, we examined how DNAmTL versus qPCR-measured TL (qPCR-TL) correlated with measures typically associated with TL. Contrary to previous findings, on almost all measures of external validity (correlations with parental TLs, southern blot TL, and age), qPCR-TL outperformed DNAmTL. The "kilobase" units of DNAm-based estimators of TL showed considerable deviations from southern blot-derived kilobase measures. Our findings suggest that DNAmTL is not a reliable index of inherited aspects of TL and underscores uncertainty about the biological meaning of DNAmTL.

Sex differences in epigenetic age in Mediterranean high longevity regions.

Sex differences in aging manifest in disparities in disease prevalence, physical health, and lifespan, where women tend to have greater longevity relative to men. However, in the Mediterranean Blue Zones of Sardinia (Italy) and Ikaria (Greece) are regions of centenarian abundance, male-female centenarian ratios are approximately one, diverging from the typical trend and making these useful regions in which to study sex differences of the oldest old. Additionally, these regions can be investigated as examples of healthy aging relative to other populations. DNA methylation (DNAm)-based predictors have been developed to assess various health biomarkers, including biological age, Pace of Aging, serum interleukin-6 (IL-6), and telomere length. Epigenetic clocks are biological age predictors whose deviation from chronological age has been indicative of relative health differences between individuals, making these useful tools for interrogating these differences in aging. We assessed sex differences between the Horvath, Hannum, GrimAge, PhenoAge, Skin and Blood, and Pace of Aging predictors from individuals in two Mediterranean Blue Zones and found that men displayed positive epigenetic age acceleration (EAA) compared to women according to all clocks, with significantly greater rates according to GrimAge (ß = 3.55; p = 1.22 × 10-12), Horvath (ß = 1.07; p = 0.00378) and the Pace of Aging (ß = 0.0344; p = 1.77 × 10-08). Other DNAm-based biomarkers findings indicated that men had lower DNAm-predicted serum IL-6 scores (ß = -0.00301, p = 2.84 × 10-12), while women displayed higher DNAm-predicted proportions of regulatory T cells than men from the Blue Zone (p = 0.0150, 95% Confidence Interval [0.00131, 0.0117], Cohen's d = 0.517). All clocks showed better correlations with chronological age in women from the Blue Zones than men, but all clocks showed large mean absolute errors (MAE >30 years) in both sexes, except for PhenoAge (MAE <5 years). Thus, despite their equal survival to older ages in these Mediterranean Blue Zones, men in these regions remain biologically older by most measured DNAm-derived metrics than women, with the exception of the IL-6 score and proportion of regulatory T cells.

The Blood DNA Methylation Clock GrimAge Is a Robust Surrogate for Airway Epithelia Aging.

One key feature of Chronic Obstructive Pulmonary Disease (COPD) is that its prevalence increases exponentially with age. DNA methylation clocks have become powerful biomarkers to detect accelerated aging in a variety of diseases and can help prognose outcomes in severe COPD. This study investigated which DNA methylation clock could best reflect airway epigenetic age when used in more accessible blood samples. Our analyses showed that out of six DNA methylation clocks investigated, DNAmGrimAge demonstrated the strongest correlation and the smallest difference between the airway epithelium and blood. Our findings suggests that blood DNAmGrimAge accurately reflects airway epigenetic age of individuals and that its elevation is highly associated with COPD.

Peroxisome proliferator-activated receptor gamma gene variants modify human airway and systemic responses to indoor dibutyl phthalate exposure.

BACKGROUND: Single nucleotide polymorphisms (SNPs) of peroxisome proliferator-activated receptor gamma (PPAR-γ; gene: PPARG) and oxidative stress genes are associated with asthma risk. However, whether such variants modulate responses to dibutyl phthalate (DBP), a common plasticizer associated with increased asthma development, remains unknown. The purpose of this study is to investigate how SNPs in PPARG and oxidative stress genes, as represented by two separate genetic risk scores, modify the impact of DBP exposure on lung function and the airway and systemic response after an inhaled allergen challenge. METHODS: We conducted a double-blinded human crossover study with sixteen allergen-sensitized participants exposed for three hours to DBP and control air on distinct occasions separated by a 4-week washout. Each exposure was followed by an allergen inhalation challenge; subsequently, lung function was measured, and blood and bronchoalveolar lavage (BAL) were collected and analyzed for cell counts and allergen-specific immunoglobulin E (IgE). Genetic risk scores for PPAR-γ (P-GRS; weighted sum of PPARG SNPs rs10865710, rs709158, and rs3856806) and oxidative stress (OS-GRS; unweighted sum of 16 SNPs across multiple genes) were developed, and their ability to modify DBP effects were assessed using linear mixed-effects models. RESULTS: P-GRS and OS-GRS modified DBP effects on allergen-specific IgE in blood at 20 h (interaction effect [95% CI]: 1.43 [1.13 to 1.80], p = 0.005) and 3 h (0.99 [0.98 to 1], p = 0.03), respectively. P-GRS also modified DBP effects on Th2 cells in blood at 3 h (- 25.2 [- 47.7 to - 2.70], p = 0.03) and 20 h (- 39.1 [- 57.9 to - 20.3], p = 0.0005), and Th2 cells in BAL at 24 h (- 4.99 [- 8.97 to - 1.01], p = 0.02). An increasing P-GRS associated with reduced DBP effect on Th2 cells. Neither GRS significantly modified DBP effects on lung function parameters. CONCLUSIONS: PPAR-γ variants modulated several airway and systemic immune responses to the ubiquitous chemical plasticizer DBP. Our results suggest that PPAR-γ variants may play a greater role than those in oxidative stress-related genes in airway allergic responses to DBP. TRIAL REGISTRATION: This study reports results from The Phthalate-Allergen Immune Response Study that was registered on ClinicalTrials.gov with identification NCT02688478.

Maternal epigenetic clocks measured during pregnancy do not predict gestational age at delivery or offspring birth outcomes: a replication study in metropolitan Cebu, Philippines.

Adverse birth outcomes, such as early gestational age and low birth weight, can have lasting effects on morbidity and mortality, with impacts that persist into adulthood. Identifying the maternal factors that contribute to adverse birth outcomes in the next generation is thus a priority. Epigenetic clocks, which have emerged as powerful tools for quantifying biological aging and various dimensions of physiological dysregulation, hold promise for clarifying relationships between maternal biology and infant health, including the maternal factors or states that predict birth outcomes. Nevertheless, studies exploring the relationship between maternal epigenetic age and birth outcomes remain few. Here, we attempt to replicate a series of analyses previously reported in a US-based sample, using a larger similarly aged sample (n = 296) of participants of a long-running study in the Philippines. New pregnancies were identified prospectively, dried blood spot samples were collected during the third trimester, and information was obtained on gestational age at delivery and offspring weight after birth. Genome-wide DNA methylation was assessed with the Infinium EPIC array. Using a suite of 15 epigenetic clocks, we only found one significant relationship: advanced age on the epigenetic clock trained on leptin predicted a significantly earlier gestational age at delivery (ß = - 0.15, p = 0.009). Of the other 29 relationships tested predicting gestational age and offspring birth weight, none were statistically significant. In this sample of Filipino women, epigenetic clocks capturing multiple dimensions of biology and health do not predict birth outcomes in offspring.

Risk-focused differences in molecular processes implicated in SARS-CoV-2 infection: corollaries in DNA methylation and gene expression.

BACKGROUND: Understanding the molecular basis of susceptibility factors to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is a global health imperative. It is well-established that males are more likely to acquire SARS-CoV-2 infection and exhibit more severe outcomes. Similarly, exposure to air pollutants and pre-existing respiratory chronic conditions, such as asthma and chronic obstructive respiratory disease (COPD) confer an increased risk to coronavirus disease 2019 (COVID-19). METHODS: We investigated molecular patterns associated with risk factors in 398 candidate genes relevant to COVID-19 biology. To accomplish this, we downloaded DNA methylation and gene expression data sets from publicly available repositories (GEO and GTEx Portal) and utilized data from an empirical controlled human exposure study conducted by our team. RESULTS: First, we observed sex-biased DNA methylation patterns in autosomal immune genes, such as NLRP2, TLE1, GPX1, and ARRB2 (FDR < 0.05, magnitude of DNA methylation difference Δß > 0.05). Second, our analysis on the X-linked genes identified sex associated DNA methylation profiles in genes, such as ACE2, CA5B, and HS6ST2 (FDR < 0.05, Δß > 0.05). These associations were observed across multiple respiratory tissues (lung, nasal epithelia, airway epithelia, and bronchoalveolar lavage) and in whole blood. Some of these genes, such as NLRP2 and CA5B, also exhibited sex-biased gene expression patterns. In addition, we found differential DNA methylation patterns by COVID-19 status for genes, such as NLRP2 and ACE2 in an exploratory analysis of an empirical data set reporting on human COVID-9 infections. Third, we identified modest DNA methylation changes in CpGs associated with PRIM2 and TATDN1 (FDR < 0.1, Δß > 0.05) in response to particle-depleted diesel exhaust in bronchoalveolar lavage. Finally, we captured a DNA methylation signature associated with COPD diagnosis in a gene involved in nicotine dependence (COMT) (FDR < 0.1, Δß > 0.05). CONCLUSION: Our findings on sex differences might be of clinical relevance given that they revealed molecular associations of sex-biased differences in COVID-19. Specifically, our results hinted at a potentially exaggerated immune response in males linked to autosomal genes, such as NLRP2. In contrast, our findings at X-linked loci such as ACE2 suggested a potentially distinct DNA methylation pattern in females that may interact with its mRNA expression and inactivation status. We also found tissue-specific DNA methylation differences in response to particulate exposure potentially capturing a nitrogen dioxide (NO2) effect-a contributor to COVID-19 susceptibility. While we identified a molecular signature associated with COPD, all COPD-affected individuals were smokers, which may either reflect an association with the disease, smoking, or may highlight a compounded effect of these two risk factors in COVID-19. Overall, our findings point towards a molecular basis of variation in susceptibility factors that may partly explain disparities in the risk for SARS-CoV-2 infection.

DNA methylation of HPA-axis genes and the onset of major depressive disorder in adolescent girls: a prospective analysis.

The stress response system is disrupted in individuals with major depressive disorder (MDD) as well as in those at elevated risk for developing MDD. We examined whether DNA methylation (DNAm) levels of CpG sites within HPA-axis genes predict the onset of MDD. Seventy-seven girls, approximately half (n = 37) of whom were at familial risk for MDD, were followed longitudinally. Saliva samples were taken in adolescence (M age = 13.06 years [SD = 1.52]) when participants had no current or past MDD diagnosis. Diagnostic interviews were administered approximately every 18 months until the first onset of MDD or early adulthood (M age of last follow-up = 19.23 years [SD = 2.69]). We quantified DNAm in saliva samples using the Illumina EPIC chip and examined CpG sites within six key HPA-axis genes (NR3C1, NR3C2, CRH, CRHR1, CRHR2, FKBP5) alongside 59 genotypes for tagging SNPs capturing cis genetic variability. DNAm levels within CpG sites in NR3C1, CRH, CRHR1, and CRHR2 were associated with risk for MDD across adolescence and young adulthood. To rule out the possibility that findings were merely due to the contribution of genetic variability, we re-analyzed the data controlling for cis genetic variation within these candidate genes. Importantly, methylation levels in these CpG sites continued to significantly predict the onset of MDD, suggesting that variation in the epigenome, independent of proximal genetic variants, prospectively predicts the onset of MDD. These findings suggest that variation in the HPA axis at the level of the methylome may predict the development of MDD.

Integrated analysis of environmental and genetic influences on cord blood DNA methylation in new-borns.

Epigenetic processes, including DNA methylation (DNAm), are among the mechanisms allowing integration of genetic and environmental factors to shape cellular function. While many studies have investigated either environmental or genetic contributions to DNAm, few have assessed their integrated effects. Here we examine the relative contributions of prenatal environmental factors and genotype on DNA methylation in neonatal blood at variably methylated regions (VMRs) in 4 independent cohorts (overall n = 2365). We use Akaike's information criterion to test which factors best explain variability of methylation in the cohort-specific VMRs: several prenatal environmental factors (E), genotypes in cis (G), or their additive (G + E) or interaction (GxE) effects. Genetic and environmental factors in combination best explain DNAm at the majority of VMRs. The CpGs best explained by either G, G + E or GxE are functionally distinct. The enrichment of genetic variants from GxE models in GWAS for complex disorders supports their importance for disease risk.

A biologically-informed polygenic score identifies endophenotypes and clinical conditions associated with the insulin receptor function on specific brain regions.

BACKGROUND: Activation of brain insulin receptors modulates reward sensitivity, inhibitory control and memory. Variations in the functioning of this mechanism likely associate with individual differences in the risk for related mental disorders (attention deficit hyperactivity disorder or ADHD, addiction, dementia), in agreement with the high co-morbidity between insulin resistance and psychopathology. These neurobiological mechanisms can be explored using genetic studies. We propose a novel, biologically informed genetic score reflecting the mesocorticolimbic and hippocampal insulin receptor-related gene networks, and investigate if it predicts endophenotypes (impulsivity, cognitive ability) in community samples of children, and psychopathology (addiction, dementia) in adults. METHODS: Lists of genes co-expressed with the insulin receptor in the mesocorticolimbic system or hippocampus were created. SNPs from these genes (post-clumping) were compiled in a polygenic score using the association betas described in a conventional GWAS (ADHD in the mesocorticolimbic score and Alzheimer in the hippocampal score). Across multiple samples (n = 4502), the biologically informed, mesocorticolimbic or hippocampal specific insulin receptor polygenic scores were calculated, and their ability to predict impulsivity, risk for addiction, cognitive performance and presence of Alzheimer's disease was investigated. FINDINGS: The biologically-informed ePRS-IR score showed better prediction of child impulsivity and cognitive performance, as well as risk for addiction and Alzheimer's disease in comparison to conventional polygenic scores for ADHD, addiction and dementia. INTERPRETATION: This novel, biologically-informed approach enables the use of genomic datasets to probe relevant biological processes involved in neural function and disorders. FUND: Toxic Stress Research network of the JPB Foundation, Jacobs Foundation (Switzerland), Sackler Foundation.